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Abstract 

We cross-correlate the Cosmic Microwave Background (CMB) temperature anisotropics observed 
by the Wilkinson Microwave Anisotropy Probe (WMAP) with the projected distribution of ex- 
tended sources in the Two Micron All Sky Survey (2MASS). By modelling the theoretical expecta- 
tion for this signal, we extract the signatures of dark energy (Integrated Sachs- Wolfe effect;ISW), 
hot gas (thermal Sunyaev-Zeldovich effect ;thermal SZ), and microwave point sources in the cross- 
correlation. Our strongest signal is the thermal SZ, at the 3.1 — 3.7a level, which is consistent with 
the theoretical prediction based on observations of X-ray clusters. We also see the ISW signal at 
the 2.5a level, which is consistent with the expected value for the concordance ACDM cosmology, 
and is an independent signature of the presence of dark energy in the universe. Finally, we see the 
signature of microwave point sources at the 2.7a level. 
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I. INTRODUCTION 

The recently released WMAP P| results constrain our cosmology with an unprecedented 
accuracy. Most of these constraints come from the linear fossils of the early universe which 
have been preserved in the temperature anisotropics of the CMB. These are the ones that can 
be easily understood and dealt with, within the framework of linear perturbation theory. 
However, there are also imprints of the late universe which could be seen in the WMAP 
results. Most notably, the measurement of the optical depth to the surface of last scattering, 
r ~ 0.17, which implied an early reionization of the universe, was the biggest surprise. There 
is also the strangely small amplitude of the large-angle CMB anisotropics which remains 
unexplained 13]. 

Can we extract more about the late universe from WMAP? Various secondary effects 
have been studied in the literature (see e.g. 23[ which lists a few). The main secondary 
anisotropy at large angles is the so-called Integrated Sachs- Wolfe (ISW) effect j^sl, which is 
a signature of the decay of the gravitational potential at large scales. This could be either a 
result of spatial curvature, or presence of a component with negative pressure, the so-called 
dark energy, in the universe j40|]. Since WMAP has constrained the deviation from flatness 
to less than 4%, the ISW effect may be interpreted as a signature of dark energy. At smaller 
angles, the dominant source of secondary anisotropy is the thermal Sunyaev-Zeldovich (SZ) 
effect |5^], which is due to scattering of CMB photons by hot gas in the universe. 

However, none of these effects can make a significant contribution to the CMB power 
spectrum below I ~ 1000, and thus they are undetectable by WMAP alone. One pos- 
sible avenue is cross-correlating CMB anisotropics with a tracer of the density in the late 



universe! 



m 



41|. This was first done by 8 who cross-correlated the COBE/DMR map 
with the NRAO VIA Sky Survey (NVSS)|lQ|. After the release of WMAP, different groups 
cross-correlated the WMAP temperature maps with various tracers of the low-redshift uni- 
verse. This was first done with the ROSAT X-ray map in [14], where a non-detection of the 
thermal SZ effect puts a constraint on the hot gas content of the universe. claimed a 
2-5(7 detection of an SZ signal by filtering WMAP maps via templates made using known 
X-ray cluster catalogs. |.38i] looked at the cross-correlation with the NVSS radio galaxy sur- 
vey, while 0] repeated the exercise for NVSS, as well as the HEAO-1 hard X-ray background 
survey, both of which trace the universe around redshift of ~ 1. Both groups found their 
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result to be consistent with the expected ISW signal for the WMAP concordance cosmology 
P], i.e. a flat ACDM universe with Qm — 0.3, at the 2 — 3cr level. Their result is consistent 
with the ACDM paradigm which puts most of the energy of the universe in dark energy 4^ . 

3J] which traces 



More recently, cross-correlated WMAP with the APM galaxy survey 
the galaxy distribution at 2; ~ 0.15. This led to an apparent detection of both the thermal 
SZ and ISW signals. However, the fact that they use a jack-knife covariance matrix to 
estimate the strength of their signal, while their jack-knife errors are significantly smaller 
than those obtained by Monte-Carlo realizations of the CMB sky (compare their Figure 
2 and Figure 3) weakens the significance of their claim. Indeed, as we argue below (see 
III.C), using Monte-Carlo realizations of the CMB sky is the only reliable way to estimate 
a covariance matrix if the survey does not cover the whole sky. 

[3^ cross-correlates the highest frequency band (W-band) of the WMAP with the ACO 
cluster survey as well as the galaxy groups and clusters in the APM galaxy survey. 
They claim a 2.6a detection of temperature decrement on angles less than 0.5°, which they 
associate with the thermal SZ effect. However they only consider the Poisson noise in their 
cluster distribution as their source of error. This underestimates the error due to large spatial 
correlations (or cosmic variance) in the cluster distribution (see III.C). also studies the 
cross-correlation of the W-band with the NVSS radio sources below a degree and claims 
a positive correlation at the scale of the W-band resolution. This may imply a possible 
contamination of the ISW signal detection in l9| and 138|. However the achromatic nature 

Q 

of this correlation makes this unlikely LD^. 

n n 

Finally, |4a| and [12| repeated the cross-correlation analysis with the 3400 and 2000 square 
degrees, respectively, of the Sloan Digital Sky Survey |52]. Both groups claim detection of a 
positive signal, but they both suffer from the inconsistency of their jack-knife and Monte- 
Carlo errors. 

The 2MASS Extended Source Catalog (XSC)|26^ is a full sky, near infrared survey of 
galaxies whose median redshift is around z ~ 0.1. The survey has reliable and uniform 
photometry of about 1 million galaxies, and is complete, at the 90% level for K-magnitudes 
brighter than 14, over ~ 70% of the sky. The large area coverage and number of galaxies 
makes the 2MASS XSC a good tracer of the ISW and SZ signals in the cross-correlation 
with the CMB. 

In this paper, we study the cross-correlation of the WMAP Q,V and W bands 
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with four different K-magnitude bins of the 2MASS Extended Source Catalog, and 
fit it with a three component model which includes the ISW, thermal SZ effects and 
microwave sources. We compare our findings with the theoretical expectations from 
the WMAP+CBI+ACBAR+2dF+Ly-a best fit cosmological model (WMAP concordance 

n 

model from here on; see Table 3 in [5|]), which is a flat universe with, = 0.27, Qb = 
0.044, h = 0.71, and as = 0.84. We also assume their values of Ug = 0.93, and 
dug/dlnk = —0.031 for the spectral index and its running at k = 0.05 Mpc j49| . 

We briefly describe the relevant secondary anisotropics of the CMB in Sec. II. Sec. Ill 
describes the properties of the cross-correlation of two random fields, projected on the sky. 
Sec. IV summarizes the relevant information on the WMAP temperature maps and the 
2MASS Extended Source Catalog. Sec. V describes our results and possible systematics, 
and Sec. VI concludes the paper. 



II. WHAT ARE THE SECONDARY ANISOTROPIES? 

The dominant nature of the Cosmic Microwave Background (CMB) fluctuations, at angles 
larger than ~ 0.1 degree, is primordial, which makes CMB a snapshot of the universe at 
radiation-matter decoupling, around redshift of ~ 1000. However, a small part of these 
fluctuations can be generated as the photons travel through the low redshift universe. These 
are the so-called secondary anisotropics. In this section, we go through the three effects which 
should dominate the WMAP/2MASS cross-correlation. 



A. Integrated Sachs- Wolfe effect 



The first one is the Integrated Sachs- Wolfe (ISW) effect 43| which is caused by the time 
variation in the cosmic gravitational potential, $. For a flat universe, the anisotropy due to 
the ISW effect is an integral over the conformal time r] 

(5isw(n) = 2 J <i>'[{r]o-r])h,r]] dr], (1) 

where $' = d^/dr], and n is unit vector in the line of sight. The linear metric is assumed 
to be 

ds'^ = a\ri){[l + 2<l>(x, ri)]dri'^ - [1 - 2<I)(x, r/)]dx ■ dx}, (2) 
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and ?7o is the conformal time at the present. 

In a flat universe, the gravitational potential $ is constant for a fixed equation of state 
and therefore observation of an ISW effect is an indicator of a change in the equation of 
state of the universe. Assuming that this change is due to an extra component in the matter 
content of the universe, the so-called dark energy, this component should have a negative 
pressure to become important at late timesQ. Therefore, observation of an ISW effect in 
a flat universe is a signature of dark energy. 

The ISW effect is observed at large angular scales because most of the power in the 
fluctuations of $ is at large scales. Additionally, the fluctuations at small angles tend to 
cancel out due to the integration over the line of sight. 

B. Thermal Sunyaev-Zeldovich effect 

The other significant source of secondary anisotropies is the so-called thermal Sunyaev- 
Zeldovich (SZ) effect js^, which is caused by the scattering of CMB photons off the hot 
electrons of the intra-cluster medium. This secondary anisotropy is frequency dependent, 
i.e. it cannot be associated with a single change in temperature. If we define a frequency 
dependent T(z/) so that Ib[v'-, ^{i')] = where /(z/) is the CMB intensity and IbIi^', T] is 
the black-body spectrum at temperature T, the SZ anisotropy takes the form 



and Pe is the electron pressure. Assuming a linear pressure bias with respect to the matter 
over density 6m- 





where 




(4) 



6pe 




(5) 



Pe 



Eq.(3) can be written as 





where 



T 





(7) 
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and Te and He are the average temperature and the comoving density of (all) electrons, 
respectively. In Appendix A, we make an analytic estimate for Tg, based on the mass 
function and mass-temperature relation of galaxy clusters. 



C. Microwave Sources 



Although technically they are not secondary anisotropies, microwave sources may con- 
tribute to the cross-correlation signal, as they are potentially observable by both WMAP 
and 2MASS. For simplicity, we associate an average microwave luminosity with all 2MASS 
sources. We can relax this assumption by taking this luminosity to be a free parameter for 
each magnitude bin, and/or removing the clustering of the point sources. As we discuss in 
Sec. V.C, neither of these change our results significantly. 

For the microwave spectrum in different WMAP frequencies, we try both a steeply falling 
antenna temperature oc l/z/^~^ (consistent with WMAP point sources ji^) and a Milky Way 
^pe spectrum which we obtain from the WMAP observations of the Galactic foreground 

In Appendix B, assuming an exponential surface emissivity with a scale length of 5 kpc 
for the Galactic disk and a small disk thickness, we use the Galactic latitude dependence of 
the WMAP temperature to determine the luminosity of the Milky Way (Eq.B6) in different 
WMAP bands: 

L*Q = 1.7 X 10^^ ergs-\ 
L*y = 3.0 X 10^^ ergs~\ 
and L^j, = 1.0 X 10^^ ergs-\ (56) 
These values are within 50% of the observed WMAP luminosity of the Andromeda 



galaxy (see Appendix B) [55[. In V.C, we compare the observed average luminosity of the 
2MASS sources to these numbers (see Table II). 

The contribution to the CMB anisotropy due to Point Sources (see Eq.B2) is given by 

where Ax is the effective bandwidth of the WMAP band|3], ndr) is the average comov- 
ing number density of the survey galaxies, di is luminosity distance, and 6g is the galaxy 
overdensity. 



III. THE CROSS-CORRELATION POWER SPECTRUM 



A. The Expected Signal 

We first develop the theoretical expectation value of the cross-correlation of two random 
fields, projected on the sky. Let us consider two random fields A(x) and -B(x) with their 
Fourier transforms defined as 

^k = y d^x e-* "^(x), andSk = J d^x e-^'' "S(x). (9) 

The cross-correlation power spectrum, PAsik) is defined by 

{A^M = (27r)=^53(ki-k2)PAB(A:i). (10) 

The projections of A and B on the sky are defined using Fa and Fb projection kernels 

i(n) ^ j FA{r)A{rh), and S(n) ^ j d'^ FB{r)B{rn). (11) 

For the secondary temperature anisotropics, these kernels were given in Eqs.1,6 and 8. 
For the projected galaxy overdensity, this kernel is 

F^ir) = r,1'';/\ (12) 
J ar r nc[r) 

For our treatment, we assume a constant galaxy bias, bg, which relates the galaxy fluctua- 
tions, 6g, to the overall matter density fluctuations 6m, up to a shot noise 6p 

Sg ^ bgSm + Sp. (13) 

In this work, we constrain the galaxy bias, bg, by comparing the auto-correlation of the 
galaxies with the expected matter auto-correlation in our cosmological model. Our bias, 
therefore, is model dependent. 

Now, expanding A and B in terms of spherical harmonics, the cross-power spectrum, 
Cab{^) is defined as 

CAB{i) = {AemBU 
^ f dndr2FA{ri)FB{r2)x 
/0PAi.(^)(47r)2j,(A;ri)j.(A:r2)r,™(k)F4(k) 
^ J dndr2FA(n)FB(r2) J ^Mkn)Ukr2)PAB{k), 

(14) 



where jVs are the spherical Bessel functions of rank i and Y^m's are the spherical harmonics. 

To proceed further, we use the small angle (large I) approximation for the spherical Bessel 
functions 



vr ,„ I , 2^ 



je{x) = y ^7^[^Di-(^ + 2 " + ^^^^ 

which yields 

Cab{() = j ^FAir)FB{r)PAB [^^^) " [1 + Oir% (16) 

This is the so called Limber equation 3^. As we do not use the quadrupole due to its 
large Galactic contamination, the smallest value of £ that we use is 3. Direct integration 
of Eq.(15) (for the ISW signal which is dominant for low €s, see Figure 7) shows that the 
Limber equation overestimates the cross-power by less than 2-3% at £ = 3, which is negligible 
compared to the minimum cosmic variance error (about 40%, see III.B) at this multipole. 
Therefore, the Limber equation is an accurate estimator of the theoretical power spectrum. 
Now we can substitute the results of Sec. II (Eqs.l, 6, 8 and 12) into Eq.(16) which yields 

CgAxJ) = , , \ . . / dr n,(r){2P^,^ (^.±lll\ 

47^'h^c'smh\x/2)bgL{x)] ^ / f+l/2 V ^^^^ 



where P{k) is the matter power spectrum, z is the redshift, and x is defined in Eq.(4). The 
terms in Eq.(17) are the ISW, SZ and Point Source contributions respectively. Since the 
ISW effect is only important at large scales, the cross-power of the gravitational potential 
derivative with matter fluctuations can be expressed in terms of the matter power spectrum, 
using the Poisson equation and linear perturbation theory, and thus we end up with 

^2 , 4:7i^h^c^smh'^{x/2)bgL{x) f£ + l/2 



[xk bT^^J^ Ax{l + zy \ r J 
where g is the linear growth factor of the gravitational potential, $, and g' is its derivative 
with respect to the conformal time. We will fit this model to our data in Sec. V, allowing 
a free normalization for each term. 

Finally, we write the theoretical expectation for the projected galaxy auto-power, Cgg, 
which we use to find the galaxy bias. Combining Eqs.12,13 and 16, we arrive at 

Jdr nl{r)[bl-P('-±^] + ^ ■ n-\r)] 

CggW = ^ ' \ , (19) 

[J dr nc(r)] 
8 



where the term is the power spectrum of the Poisson noise, Sp, while the extra free 
parameter, 7, is introduced to include the possible corrections to the Poisson noise due to 
the finite pixel size. In the absence of such corrections 7 = 1. In Sec. V, we seek the values 
of bg and 7 that best fit our observed auto-power for each galaxy sub-sample. 

To include the effects of non-linearities in the galaxy power spectrum, we use the Peacock 



& Dodds fitting formula 



for the non-linear matter power spectrum, P{k). 



B. Theoretical errors: cosmic variance vs. shot noise 

To estimate the expected theoretical error, again for simplicity, we restrict the calculation 
to the small angle limit. In this limit, the cross-correlation function can be approximated 

by 

CABii) ^ ^{A,^ b;J, (20) 

where AQ is the common solid angle of the patch of the sky covered by observations of both 
A and 50. 

Assuming gaussianity, the standard deviation in Cab, for a single harmonic mode, is 
given by 

ACUi) = {CUi)) - {CABii))' 

= AQ'^[{Aim B^^){Aim B}^) + {Aim AiJ){Bim ^L)] 

= Cl^W + CAAii)CBBii). (21) 

The number of modes available between i and £ + 1, in the patch AQ, is 

ATV^P^y)^" (22) 

4:71 

and so the standard deviation of Cab, averaged over all these modes is 

ACl^(^) ^ + CAAii)CBBm (23) 

In fact, since the main part of CMB fluctuations is of primordial origin, the first term 
in brackets is negligible for the cross- correlation error, so the error in the cross-correlation 
function, as one may expect, depends on the individual auto-correlations. 
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We can use the CMBFAST code j47| to calculate the auto-correlation of the CMB tem- 
perature fluctuations. Also, the theoretical expectation for the auto-power of the projected 
galaxy distribution is given by Eq.(19). 

The galaxy/CMB auto-power spectra are dominated by Poisson(shot) noise/detector 
noise at large £'s. Therefore, the measurement of the thermal SZ signal, which becomes 
important at large £'s, is limited by the number of observed galaxies, as well as the resolu- 
tion of the CMB detector (the angle at which signal-to-noise ratio for the CMB measurement 
is of order unity). On the other hand, for the small C. portion of the cross-correlation which 
is relevant for the ISW signal, the error is set by the matter and CMB power spectra and 
thus, is only limited by cosmic variance. The only way to reduce this error is by observing a 
larger volume of the universe in the redshift range < 2; < 1, where dark energy dominates. 



C. A Note On the Covariance Matrices 



We saw in Section III.C that the errors in cross-correlations could be expressed in terms 
of the theoretical auto-correlation. However, this is not the whole story. 

We have a remarkable understanding of the auto-power spectrum of the CMB. However, 
if one tries to use the frequency information to, say subtract out the microwave sources, the 
simple temperature auto-power does not give the cross-frequency terms in the covariance 
matrix. In fact, in the absence of a good model, the only way to constrain these terms is 
by using the cross-correlation of the bands themselves. Of course, this method is limited by 
cosmic variance and hence does not give an accurate result at low multipoles. To solve this 
problem, we use the WMAP concordance model CMB auto-power for ^ < 13. Since there is 
no frequency-dependent signal at low £'s, we only use the W-band information, which has 
the lowest Galactic contamination jo^, for our first 4 £-bins which cover 3 < £ < 13 (see the 
end of Sec. III.D for more on our £-space binning). 

There is a similar situation for the contaminants of the 2MASS galaxy catalog. Systematic 
errors in galaxy counts, due to stellar contamination or variable Galactic extinction, as 
well as observational calibration errors, may introduce additional anisotropies in the galaxy 
distribution which are not easy to model. Again, the easiest way to include these systematics 
in the error is by using the auto-correlation of the observed galaxy distribution, which is 
inaccurate for low multipoles, due to cosmic variance. Unfortunately, this is also where we 
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expect to see possible Galactic contamination or observational systematics. With this in 
mind, we try to avoid this problem by excluding the quadrupole, C(2), from our analysis. 



At this point, we should point out a misconception about the nature of Monte-Car 



16 



o vs. 

[21. 



jack-knife error estimates in some previous cross-correlation analyses, specifically 
Many authors have used Gaussian Monte-Carlo realizations of the CMB sky to estimate 
the covariance matrix of their real-space cross-correlation functions 

S Q, H Q- The 

justification for this method is that, since the first term in Eq.(23) is much smaller than 
the second term, the error in cross-correlation for any random realization of the maps is 
almost the same as the true error, and the covariance of the cross-correlation, obtained from 
many random Gaussian realizations is an excellent estimator of the covariance matrix. We 
may also obtain error estimates based on random realizations of one of the maps, as long 
as the observed auto-power is a good approximation of the true auto-power, i.e. the cosmic 
variance is low, which should be the case for angles smaller than 20 degrees {C, > 10). Of 
course, at larger angles, as we mentioned above, one is eventually limited by the systematics 
of the galaxy survey and, unless they are understood well enough, since theoretical error 
estimate is not possible, there will be no better alternative rather than Monte-Carlo error 



estimates. In fact, contrary to 



ss-correlation makes 

Monte-Carlo errors a slight underestimate (see Eq.23). On the other hand, there is no 
rigorous justification for the validity of jack-knife covariance matrices, and the fact that 



Hk-knife errors could be smaller than the Monte-Carlo errors by up to a factor of three 

As we argue below (see Sec. III.D), since we do our analyses in harmonic space and 
use most of the sky, our £-bins are nearly independent and performing Monte-Carlo's is 
not necessary. Our covariance matrix is nearly diagonal in £-space and its elements can be 
obtained analytically, using Eq.(23). 



D. Angular Cross-Power Estimator 

The WMAP temperature maps are set in HEALPix format which is an equal area, 
iso-latitude pixellization of the sky. As a component of the HEALPix package, the EFT 
based subroutine 'map2alm' computes the harmonic transform of any function on the whole 
sky. However, as we describe in the next section, in order to avoid contamination by Galactic 
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foreground emission in WMAP temperature maps, and contamination by stars and Galactic 
extinction in the 2MASS survey, we have to mask out ~ 15% of the CMB and ~ 30% of the 
2MASS sky. Therefore, we cannot obtain the exact values of the multipoles, and should 
use an estimator. 

We use a quadratic estimator which is based on the assumption that our masks, W{n), 
are independent of the data that we try to extract (see j3| for a review of different es- 
timators). The real-space cross-correlation of the masked fields ^(n) = WA{ii)A{fi) and 
-B(n) = WB{ii)B{h) on the sphere is given by 

{A{h)B{m)){A{h)B{m)WA{n)WB{ni)) 

= {A{n)B{ih)){WA{n)WB{ni)), (24) 

where, in the last step, we used the independence of data and masks, and averaged over 
all pairs of the same separation. Assuming that (ll^A(n)lFB(m)) does not vanish for any 
separation (which will be true if the masked out area is not very large), we can invert this 
equation and take the Legendre transform to obtain the un-masked multipoles 



CAsi^) = Y.F,eCAB{n, where 



1, f Pe(cose)Pe'(cose) , 



In 

team 



act this estimator is mathematically identical 



to the one used by the WMAP 



2l| . and, within the computational error, should give the same result. The difference 
is that we do the inversion in real-space, where it is diagonal, and then transform to harmonic 
space, while they do the inversion directly in harmonic space. Indeed, using our method, we 
reproduce the WMAP binned multipoles [2^ within 5%. However, we believe our method 
is computationally more transparent and hence more reliable. Also, the matrix inversion in 
harmonic space is unstable for a small or irregular sky coverage (although it is not relevant 
for our analyses). 

Finally, we comment on the correlation among different multipoles in £-space. Masking 
about 30% of the sky causes about 30% correlation among neighboring multipoles. We 
bin our multipoles into 13 bins that are logarithmically spaced in i (covering 3 < i < 
1000) , while excluding the quadrupole due to its large Galactic contamination in both data 
sets. The highest correlation between neighboring bins is 15% between the first and the 
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second bins (C(3) and (9C(4) + llC(5))/20). To simplify our calculations, we neglect this 
correlation, as any correction to our results will be significantly smaller than the cosmic 
variance uncertainty(see V.B), i.e. we approximate our covariance matrix as diagonal in 
£-space. 



IV. DATA 



A. WMAP foreground cleaned temperature maps 

We use the first year of the observed CMB sky by WMAP for our analysis jsj. The 
WMAP experiment observes the microwave sky in 5 frequency bands ranging from 23 to 
94 GHz. The detector resolution increases monotonically from 0.88 degree for the low- 
est frequency band to 0.22 degree for the highest frequency. Due to their low resolution 
and large Galactic contamination, the two bands with the lowest frequencies, K(23 GHz) 
and Ka(33 GHz), are mainly used for Galactic foreground subtraction and Galactic mask 
construction |6|, while the three higher frequency bands, which have the highest resolution 
and lowest foreground contamination, Q(41 GHz), V(61 GHz), and W(94 GHz), are used for 
CMB anisotropy spectrum analysis, use the Maximum Entropy Method to combine the 
frequency dependence of 5 WMAP bands with the known distribution of different Galactic 
components that trace the dominant foregrounds to obtain the foreground contamination 
in each band. This foreground map is then used to clean the Q, V and W bands for the 
angular power spectrum analysis. Similarly, we use the cleaned temperature maps of these 
three bands for our cross-correlation analysis. We also use the same sky mask that they use, 
the Kp2 mask which masks out 15% of the sky, in order to avoid any remaining Galactic 
foreground, as well as 208 identified microwave point sources. 



B. 2MASS extended source catalog 



We use galaxies from the Near-IR Two Micron All Sky Survey [2MASS;|^ as the large- 
scale structure tracer of the recent universe. Our primary data set is the public full-sky 
extended source catalog [XSC; 26]. The f^g-band isophotal magnitude, is the default 
flux indicator we use to select the external galaxies for our analysis. i^2o is the measured flux 
inside a circular isophote with surface brightness of 20 mag arcsec"^. The raw magnitudes 
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from the catalog were corrected for Galactic extinction using the IR reddening map of 
Schlegel, Finkbeiner & Davis 45[: 



i^20 ^ i^2o - Ak, (26) 

where Ak = RkE{B — V) = 0.367 x E{B — V) [57]. There are approximately 1.5 million 
extended sources with corrected K20 < 14.3 after removing known artifacts (ccjlag != 'a' 
and 'z') and using only sources from a uniform detection threshold (use_src = 1). 



1. Completeness and Contamination 

We use the standard log N-log S test to determine the completeness limit of the extended 
source catalog. The top panel of Figure Q shows the number of galaxies as a function of K2o- 
The log number counts can be approximated by a power-law: 

^ oc 10'^'". (27) 
dm 

To infer the true number count-magnitude relation, we need to ensure that our catalog is 
free from contaminants since not all extended sources from the XSC are external galaxies. 
At low Galactic latitude where stellar density is high, unresolved multiple star systems are 
often confused as extended sources. For the purpose of fitting for the power-law slope k,, 
we use only sources with |6| > 30°. Using ~ 250,000 galaxies in the magnitude range 
.3.2 < m < 13.7 (where the reliability has been determined to be 99% by Huchra & Mader 



2J|), we fitted a number count slope k = 0.676 ± 0.005. 

While the XSC is unreliable at low Galactic latitudes, the \b\ > 30° cut is too permissive 
and would throw away a large area of the sky that can be used for analysis. In principle, 
we could use the stellar density Ustar from the 2MASS Point-Source Catalog(PSC) to set a 
;hreshold for excluding region of high stellar density. However, since it has been shown by 
3 and that unresolved extended galaxies are found in the PSC(up to 2% of all point 



sources with K20 ~ 14), a mask derived from the observed stellar density would preferentially 
exclude regions of high galaxy density. 

We use the extinction map of [43] to exclude regions of the sky where the XSC is un- 
reliable. Figure El shows the average number of galaxies per HEALPix pixel of 0.83 deg^ 
{Nside = 64), as a function of Galactic extinction for the four magnitude ranges used in 
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FIG. 1: (Top panel) The histogram is the observed K20 number-magnitude relation for galaxies 
in regions with Ak < 0.05. The solid line is the model counts inferred using data from \b\ > 30° in 
the magnitude range 13.2 < K20 < 13.7 where the extended source catalog(XSC) is most reliable. 
(Bottom panel) The square points gives the completeness as inferred from the difference between 
the observed and model counts. The solid curve is a fit to a parametric model that estimates both 
the incompleteness and contamination rate in a consistent manner. The dotted curve is a similar 
fit using data with a less stringent Ak < 0.2. The dashed curve is from |6| > 30°, which serves 
roughly as the completeness upper-bound for the XSC. The vertical line at K20 = 13.85 gives a 
completeness at 98% for data with Ak < 0.05 used in our analysis. 
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FIG. 2: Average number of galaxies per 0.83 deg^ pixel (HEALpix Nside = 64) as a function of 
extinction. For bright galaxies {K20 < 13.5), the galaxy density is constant up to extinction value 
~ 0.25. For 13.5 < K20 < 14.0, the density drops off at Ax ~ 0.65. We use only regions with 
Ak < 0.05 {dashed vertical line) for our analysis. Errors are estimated using jack-knife resampling. 



our analysis. For bright galaxies, e.g. K20 < 13.5, the Galactic density is constant on de- 
gree scales. For the faintest magnitude bin, the number density drops off at large for 
Ak beyond ~ 0.065. We thus choose A^ < 0.05 j^]. This stringent threshold excludes 
~ 99% of all regions with Ugtar > 5000 deg^^. Moreover, it improves the global reliability of 
galaxy counts, as our flux indicator K20 for each source was corrected for Galactic extinc- 
tion, which has an uncertainty that scales with Ak itself. This cut reduces the number of 
extended sources with K20 < 14.3 to ~ 1 million, covering ~ 68.7% of the sky. For the sake 
of completeness, we also repeat our cross-correlation analysis for a less stringent mask with 
Ak < 0.1, which covers ~ 79.0% of the sky. 

Using K = 0.676 derived from regions with |6| > 30° as a model for the true underlying 
number counts, we infer the catalog completeness and contamination as a function of ap- 
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parent magnitude for the extinction cropped sky. We deduce the intercept of the hnear log 
counts - magnitude model by scaling the observed number counts from the |6| > 30 region 
to the larger masked sky at the bright magnitude range 12.5 < K20 < 13.0. Essentially, 
we assumed the two number count distributions are identical at those magnitudes. The 
observed fractional deviation from Eq. (j27p 

I(m) - f — " - (28) 
\dm dm J I dm 

is positive at faint magnitudes indicating incompleteness but crosses zero to a constant 
negative level towards the bright-end, which we inferred as contamination to the XSC 
Plotted in the bottom panel of Figure ^ is the completeness function C{m) = 1 — I{m), 
where we parametrically fitted using 

(m — m)2 



/(m) = lo exp 



Const. (29) 



2^2 

In Figure H the term Const describes the low level of excursion beyond C{m) = 1. We 
obtained a ~ 98% completeness for K20 < 13.85 and contamination rate at 0.5% level 
for Ak < 0.05 (solid curve). As a comparison, a less stringent threshold of Ax < 0.2 
(dotted curve), the completeness is ~ 95% with contamination at 1.5%. The dashed curve 
is computed using high latitude data {\b\ > 30°), serves roughly as the completeness upper- 
bound (as a function of apparent magnitude) for the XSC. 

At a low level, contaminants in the catalog merely increase the noise of our signal with 
marginal systematic bias. The A^ < 0.05 extinction mask is close to optimal in terms of 
signal-to-noise for our cross-correlation analysis. One the other hand, catalog incompleteness 
at faint magnitudes affects our ability to infer the correct redshift distribution. We use 
galaxies up to K20 = 14.0 but weighted the redshift distribution at a given magnitude range 
(described below) by Eq. (|^. 



2. Redshift Distribution 



The redshift distribution of our sample was inferred from the Schechter j44| parameters 
fit of the K20 luminosity function from 2^. The redshift distribution, dN/dz of galaxies 
in the magnitude range mbright < m < mfaint is given by the integration of the luminosity 
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TABLE I: -7— parameters for the four magnitude bins 








P 


A 


AT 
iV 


12.0 < K20 < 12.5 


0.043 


1.825 


1.574 


49817 


12.5 < K20 < 13.0 


0.054 


1.800 


1.600 


102188 


13.0 < K20 < 13.5 


0.067 


1.765 


1.636 


211574 


13.5 < K20 < 14.0 


0.084 


1.723 


1.684 


435232 



function $(M) 

[z)dz= <^{M)dM x^{z)dz, (30) 



dz JMbiz) dz 

where dVc/dz is the Une-of-sight comoving volume element and 

Mf{z) ^ mfaint - DM{z) - k{z) (31) 
Mfe(2;) = mbright - DM{z) - k{z) (32) 

Here, DM{z) and k{z) are the distance modulus and /c-correction at redshift z. To be 
consistent with [29], we employ k{z) = — 6.01og(l + z), but the redshift distribution is 
insensitive to the exact form of the A;-correction. The Schechter parameters used were 
M* = —23.39 and a.^ = —1.09. For analytic convenience, we further model dN/dz with the 
three parameter generalized-gamma distribution: 

— — I A, p, Zo) dz oc 

i \ZoJ 



dz r(A) \Zo 

X exp 



z 



Zo 



The fit were weighted by relative counts, hence they are exact near Zo-, the mode of the 
distribution, but underestimate the true number at the high redshift tail by less than 1%. 
Table H] gives the redshift distribution parameters for the four magnitude bins used in our 
analysis. We normalize the integral of dN / dz with the total number of observed galaxies in 
the respective magnitude range, 

/■°° dN 

iVtotai(Al]) = / —{z)dzx^^l (34) 

^ dVc f 

nc{z) ~^{z) dz / dfl 

nc{r) dr dVL. (35) 
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FIG. 3: dN/dz for the four magnitude bins used in the analysis. 

The observed A^totai is consistent with the 10% uncertainty in the normahzation, cj)* = 
1.16 X 10~^ /i^ Mpc^"^ obtained by i22| in their luminosity function analysis. Eq. flH5|l gives 
the explicit relation (in the absence of clustering) between the redshift distribution, dN/ dz, 
and the comoving density, nc(r) used in Sec. III. A. 

Figure El is a plot of the redshift distribution for the four magnitude bins used in our 
analysis. For the first three bright samples, where we are complete, the parameters for 
dN/dz were derived from a direct application of Eq. 1213 For 13.5 < K20 < 14.0, the redshift 
distribution was computed by summing up magnitude slices with interval AK20 = 0.05, and 
weighted by their relative number counts. 



V. RESULTS 

For the following results, unless we mention otherwise, we use the WMAP concordance 
cosmological model. 

By comparing the angular auto-power spectrum of the galaxies in each magnitude bin 
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FIG. 4: The auto-power for our four different magnitude bins. The sohd curves show the observed 
auto-power multipoles with their estimated Gaussian errors (Eq.23), while the dashed curves are 



the projected Peacock and Dodds 



39| non-hnear power spectra with the best fit constant bias. 



The best fit Poisson noise term is subtracted out. 

with the theoretical auto-power spectrum (Eq.l9), we can obtain the bias of the 2MASS 
galaxies. In order to do this, we use a fit? assuming independent gaussian random errors 
at each £-bin. Figure 4 compares our best fit models of the auto-power (solid curves) with the 
measured auto-powers for each magnitude bin. The value of the bias for all the magnitude 
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bins is within 

6g = 1.11 ±0.02, (36) 

which confirms our constant bias assumption js^ . Our values for the Poisson correction factor 
(see Eq.l9), 7, are all within 1% of 1.02. The most significant deviation of the theoretical 
fit from the observed auto-power is about 30% at £ ~ 30 — 40. One possibility may be that 
galaxy bias is larger at large (linear) scales than at the (non-linear) small scales. In order 
to estimate the effect, we can limit analyses to the first 7 £-bins (£ < 70, scales larger than 
~ 7 — 13 Mpc). This yields the estimated bias on linear scales: 

6g,iin = 1.18 ±0.08, (37) 

The angular scale corresponding to £ = 30 — 40 is a few degrees, which is close to the 
length of the 2MASS scanning stripes (6°). The amplitude of deviation from the constant 
bias model would require systematic fluctuations of order 10% in the number counts on that 
scale. If these were due to systematic errors in the 2MASS photometric zero-point, such 
fluctuations would require a magnitude error Am ~ 0.06, which is significantly larger than 
the calibration uncertainties in 2MASS js^]. Therefore, we will use our estimated linear bias 
(Eq. 37) for the interpretation of our ISW signal, while we use the full bias estimate (Eq. 
36), which is dominated by non-linear scales, to analyze our SZ signal. 

The points in Figure 5 summarize our twelve observed cross-correlation functions (3 
WMAP bands x 4 magnitude bins), while Figure 6 shows the same data after subtracting out 
the best fit contribution due to microwave Point Sources. We fit our theoretical model (Eq. 
18) to our cross-correlation points (including only the W-band for the first 4 £-bins; see Sec. 
III.C), allowing for free global normalizations for the ISW, SZ and Point Source components. 
The curves show this model with the best fit normalizations for these components, while 
the shaded region shows the 68% uncertainty around a null hypothesis. Figure 7 shows how 
individual theoretical components depend on frequency and i for our faintest magnitude bin. 
As we mentioned in Sec. III.C, different £-bins are nearly independent. However, different 
combinations of frequency bands and magnitude bins are highly correlated and we use the 
full covariance matrix which we obtain from the data itself (see III.C) for our analysis. 

The apparent dispersion in our data points for the first 4-5 £-bins is smaller than what 
we expect from gaussian statistics (the shaded regions in Figures 5 and 6). This may be due 
to the non-gaussian nature of the systematics (observational or Galactic), which dominate 
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the error on large angles, and make the variance (Eq. 23) significantly different from the 
68% confidence region. 

Figures 5 & 6 show that our faintest magnitude bin has the smallest error. This is due 
to the fact that our faintest magnitude bin covers the largest effective comoving volume and 
number of galaxies (see Table I and Figure 3), and as a result, between 50-70% of our signal 
(depending on the component) comes from this bin. Repeating the statistical analysis for 
individual magnitude bins leads to results consistent with the combined analysis within the 
errors. 

A. Thermal SZ signal 

In a model described in Appendix A, we quantify the amplitude of the SZ signal in terms 
of (5, the coefficient of temperature-mass relation for clusters of galaxies 

r.(M) = (M2kaV)e(^5J5^)'". (A2) 

Different theoretical and observational methods place Q somewhere between 1 and 2, with 
observations preferring the higher end. This amplitude could be equivalcntly described in 
terms of Tg, the product of pressure bias and average electron temperature (see Eq.7), which 
is less model dependent. Our best fit value for the thermal SZ signal, which shows a signal 
at the 3.1(7 level, is 

Q = 1.19 ± 0.38, or = bpf^ = (1.04 ± 0.33) keV, (38) 

which is consistent with the X-ray observations of galaxy clusters. 

This result is slightly dependent on the spectrum of the microwave point sources, which 
we discuss below in V.C. If we restrict the analysis to i > 20, which is where all the SZ 
signal comes from, and our estimates of the covariance matrix is robust, our reduced is 
0.93 which is within the 68% allowed range for our 12 x 13 degrees of freedom. This implies 
that there is no observable systematic deviation from our theoretical expectation for the 
shape of the thermal SZ cross-power (or its gaussianity) . 

Repeating the analysis with the < 0.1 extinction mask (See Sec. I V.B.I) for the 
2MASS galaxies, which has a 10% larger sky coverage, increases the SZ signal shghtly to 
Q — 1.27 ± 0.35, which is a detection at the ~ 3.7a significance level. This is probably 
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FIG. 5: The cross-power for our four magnitude bins. The curves are the best fit model 
(ISW+SZ+Point Sources) for the three bands and the points show the data. The ISW/SZ com- 
ponents dominate the signal for ts below/above 20. The Point Source contribution becomes 
important for the lower frequency bands at the highest ts. The shaded region shows the 1 — a 
error centered at the null hypothesis. Note that, while different -^-bins are nearly independent, dif- 
ferent cross-powers of bands with magnitude bins are highly correlated. As shot noise dominates 
the signal for our last two i-h'ms, for clarity, we only show the first 11 £-bins, for which the errors 
for the three WMAP bands are almost the same. 
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FIG. 6: The same as Figure 5, but with the Point Source contributions subtracted from both 
theory and data. 

because the Galactic contamination close to the plane is only at large angles and does not 
contribute to the SZ signal. Therefore, as long as the Galactic contamination does not 
completely dominate the fluctuations, increasing the area only increases the SZ signal. 



B. ISW signal 

Using our estimated linear bias (Eq. 37), our fit yields an ISW signal of 
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FIG. 7: Different components of our best fit theoretical cross-power model, compared with the 
data for our faintest magnitude bin (13.5 < K < 14). The dotted(red) curves show the ISW 
component, while the short-dashed (green) and long-dashed (blue) curves are the SZ and Point 
Source components respectively. The black curves show the sum of the theoretical components, 
while the points are the observed cross-power data. 



ISW = 1.49 ± 0.61 (39) 
X concordance model prediction, 



a 2.5(7 detection of a cross- correlation. As with the previous cross-correlation analyses|3, 



171 |38[, this is consistent with the predictions of the concordance ACDM paradigm. 
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However, among the three signals that we try to constrain, the ISW signal is the most 
difficult to extract, because almost all the signal comes from i < 20, given our redshift 
distribution. For such low multipoles, there are several potential difficulties: 

1- The small Galactic contamination or observational systematics in 2MASS may dom- 
inate the fluctuations in the projected galaxy density at low multipoles and wipe out the 
signal. However, since we use the observed auto-power of 2MASS galaxies for our error 
estimates, this effect, which does contribute to the auto-power (and not to the signal), is 
included in our error. 

2- Our covariance estimator loses its accuracy as the cosmic variance becomes important 
at low multipoles (see III.C). A random error in the covariance matrix can systematically 
increase the ^ind hence decrease the estimated error of our signal. However, our reduced 

is 0.88, which is in fact on the low side (although within Icr) of the expected values for 
124 degrees of freedom (remember that we only used the W band for the flrst 4 ^-bins). 
Assuming gaussian statistics, this implies that we do not signiflcantly underestimate our 
error. 

3- Possible Galactic contamination in WMAP may correlate with Galactic contamination 
in 2MASS at low multipoles, which may lead to a fake positive signal. However, the largest 
contribution of Galactic foreground is visible in the Q-band j6j, and our low i multipoles 
have in fact a lower amplitude in Q band. Although this probably shows a large error due 
to contamination in the Q-band amplitude, the fact that this is lower than the amplitude of 
V and W bands implies that our main signal is not contaminated. Because of the reasons 
mentioned in lll.C, we only use the W-band information for i < 14. 

Using the less stringent extinction mask, < 0.1 (see IV.B.l), for the 2MASS sources 
yields a signal of ISW= 1.9 ± 1.1, which is a lower signal to noise detection at the 1.7a level. 
This is probably due to the fact that most of the ISW signal comes from angles larger than 
~ 10°, which is highly contaminated in regions close to the Galactic Plane. 

Finally, we should mention that since the ISW signal comes from small ts, while the SZ 
and point source signals come from large €s (See Figure 7), there is a small correlation (less 
than 10%) between the ISW and other signals. 
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TABLE II: Best fit point source strengths for different assumed spectra. The associated best fit 
SZ signal and are also quoted. Here, T4 stands for the antenna temperature, while Ly, defined 
in Eq.(B2), is the estimated luminosity of the Milky Way in WMAP's V-band. 



Spectrum 


Lv/Ly 


Q 




Milky Way 


16.2 ±7.8 


1.10 ±0.40 


111.2 


6Ta oc 


21.0 ±8.1 


1.19 ±0.38 


109.5 


6Ta oc 


10.9 ±4.7 


0.94 ± 0.33 


110.8 



C. Microwave Point Sources 

As described in II. C, we assume that our point sources trace the 2MASS objects and 
have either a Milky Way spectrum, a z/~^, or a z/~^ frequency dependence their antenna 
temperature (the last two are the expected synchrotron spectrum of radio sources). The 
results are shown in Table II. We see that, although all the spectra are consistent at a 2a level, 
we achieve the lowest for a z/~^ spectrum which is similar to the spectrum of point sources, 
identified by the WMAP team |6|. We should also note that since the ^-dependence of the 
SZ and Point Source signals are very similar, the two signals are correlated at a 50 — 70% 
level, which is shown in Figure 8. Using a less stringent extinction mask {Ak < 0.1, see 
IV.B.l) increases the detection level of microwave sources by about 10%. 

To relax our assumption for the redshift distribution of Point Sources (which we assume to 
be the same as 2MASS sources at each magnitude bin; see II. C), we can also allow different 
magnitude bins to have different Ly's and treat each as a free parameter. It turns out that 
this does not affect either our ISW or SZ signals, or their significance, while Ly's for each 
magnitude bin is consistent with the values in Table II, within the errors. As the Point 
Source signal is dominated by Poisson noise at large ^'s, removing the assumed clustering 
among the Point Sources (see II. C), does not affect the SZ or ISW signals either. 

VI. CONCLUSIONS 

We obtain the cross-power spectrum of the three highest frequency bands of the WMAP 
cosmic microwave background survey with the 2MASS Extended Source Catalog of near 
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FIG. 8: 1 and 2 — cr likelihood regions of our SZ+Point Source signals, for a STa oc i/^^ spectrum 
(see Table II). Q (defined in Eq.A2) is the coefficient of the mass-temperature relation for galaxy 
clusters, while Tg (defined in Eq.7), is the product of gas pressure bias and the average electron 
temperature. Ly is the average WMAP V-band luminosity of the 2MASS sources, while Ly 
(defined in Eq.B6) is the same number, estimated for the Milky Way. The large correlation of Q 
(SZ signal) and Ly (Point Source signal) is due to the similar ^-dependence of the two signals (see 
Figure 7) . Note that the conversion between Tg and Q depends on the assumed cosmological model 
(see Appendix A). 
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infrared galaxies. We detect an ISW signal at the ~ 2.5cr level, which confirms the presence 
of a dark energy, at a level consistent with the WMAP concordance cosmology. We also find 
evidence for an anti-correlation at small angles (large £'s), which we attribute to thermal 
SZ. The amplitude is at 3.1 — 3.7cr level and is consistent with the X-ray observations of 
galaxy clusters. Finally, we see a signal for microwave Point Sources at the 2.6ct level. 

We've seen that the completeness limit of the extended source catalog is between 13.5 
and 14 in K. However, matches with SDSS show that there are many unresolved sources 
in the 2MASS Point Source Catalog (PSC) that are in fact galaxies. If we can select out 
galaxies in the PSC, perhaps by their distinctive colors, we should be able to push the sample 
at least half a magnitude fainter than we have done here, probing higher redshifts with a 
substantially larger sample. 

Future wide-angle surveys of galaxies should be particularly valuable for cross-correlation 
with the WMAP data, especially as the latter gains signal-to-noise ratio in further data 
releases. The Pan-STARRS project [27^ for example, should yield a multi-color galaxy 
catalog to 25th mag or even fainter over 20,000 square degrees or more of the sky well before 
the end of the decade; it will more directly probe the redshift range in which the SZ and 
ISW kernels peak, and therefore should be particularly valuable for cross-correlating with 
WMAP and other CMB experiments. 
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could be the case if the CMB power is suppressed on large angles. This wou 



significance of an ISW detection at such scales. This effect is elaborated in 
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d increase the 



APPENDIX A: SEMI- ANALYTICAL ESTIMATE OF SZ SIGNAL 

In order to find Te (defined in Eq.7) we need an expression for the dependence of the 
electron pressure overdensity on the matter overdensity. As the shock-heated gas in clusters 
of galaxies has keV scale temperatures and constitutes about 5 — 10% of the baryonic mass 
of the universe, its contribution to the average pressure of the universe is significantly higher 
than the photo- ionized inter-galactic medium (at temperatures of a few eV). Thus, the 
average electron pressure in a large region of space with average density p(l + Sm) is given 
by 

= ^ [ dM-M-n{M-,p)[kB%{M-,p)][b{M) + ^^^]5rn, (Al) 
p J dlogp 

where n{M; p) and Te{M; p) are the mass function and temperature-mass relation of galaxy 
clusters respectively. Also, h{M) = ^^°gi^^'^^ is the bias factor for haloes of virial mass M 
(= M200; niass within the sphere with the overdensity of 200 relative to the critical density). 
For our analysis, we use the Sheth & Tormen analytic form for n(M) and b{M), which 
is optimized to fit numerical N-body simulations. 

We can use theoretical works on the cluster mass-temperature relation (which assume 
equipartition among thermal and kinetic energies of different components in the intra-cluster 
medium) to find Te(M), (e.g. Q) 



{0.32Q){2TTGHMy 



TJM) = (6.62 keV)g , , 
while 1 < g < 2 (A2) 
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for massive clusters, where H = lOO/i kms~^/ Mpc is the (local) Hubble constant. Although 
there is controversy on the value of the normalization Q (see e.g. 32] and references there 
in), Pi argue that, as long as there are no significant ongoing astrophysical feed-back or 
cooling (i.e., as long the evolution is adiabatic), the dependence on H and M should be the 
same. Combining this with the local comoving continuity equation 

3iH + 6H) = - 6m, (A3) 
P 



yields 

dlogTe 2d\ogH 2D 



(A4) 



dlogp 3 dlogp 9DH^ 
where D is the linear growth factor. 

One may think is that observations may be the most reliable way of constraining Q in 
Eq.(A2). However, almost all the observational signatures of the hot gas in the intra-cluster 
medium come from the X-ray observations which systematically choose the regions with 
high gas density. With this in mind, we should mention that while observations prefer a 
value of Q close to 1.7, numerical simulations and analytic estimates prefer values closer 
to 1.2^. For our analysis, we treat Q as a free parameter which we constrain using our 
cross-correlation data (see Sec. V.A). 

Putting all the pieces together, we end up with the following expression for Tg 

fe = (0.32 Q)(2'KGHf/^ 
xjdvfsT[v]M^"[hsT{v)-^l (A5) 
KM) = 

where a{M) is the variance of linear mass overdensity within a sphere that contains mass M 



of unperturbed density, whil e 5r. ~ 1.68 is the spherical top- hot linear growth thresho 
fsT and hsT are defined in For the WMAP concordance cosmological modeljs], this 

integral can be evaluated to give 

= fepfe = (0.88 Q) keV. (A6) 

The above simple treatment of the SZ signal fails at scales comparable to the minimum 
distance between clusters, where the average gas pressure does not follow the average matter 
density 0, S|, which leads to a scale-dependent pressure bias. Moreover, efficient galaxy 
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formation removes the hot gas from the intra-cluster medium, which causes Eq.Al to over- 
estimate the SZ signal. As this paper mainly focuses on the observational aspects of our 
detection, we delay addressing these issues into a further publication jsil. The preliminary 
results seem to be consistent with the above simple treatment at the 20% level. 



APPENDIX B: MICROWAVE LUMINOSITIES OF THE ANDROMEDA 
GALAXY AND THE MILKY WAY 

First we derive how much the flux received by a microwave source at distance di and 
observed solid angle 5VL affects the observed CMB temperature. The apparent change in the 
black-body temperature is obtained by 



5n -5 



Audi' 



(Bl) 



_exp[hu/{kBT^^,^)] - 1_ 
The left hand side of Eq.(Bl) is the change in the Planck intensity, where u and Au are 
the detector frequency and band width respectively. The right hand side is the observed 
Microwave flux. Defining x as the frequency in units of ksT^^j^/h, Eq.(Bl) yields 

6T _ An^h^c^ sinh2(x/2) L 

To obtain the microwave luminosity of Milky Way, we assume an optically and geomet- 
rically thin disk, with a microwave emissivity, e, which is constant across its thickness and 
falls as eoexp(— r/ro) with the distance, r, from its center. The disk thickness is 2H ^ r, 
while we assume tq ^ 5kpc, our distance from the Galactic center is r ~ 8.5 kpc, and our 
vertical distance from the center of the disk is z. Integrating Eq.(B2) over the disk thickness 
leads to the cosecant law for the Galactic emission 
5T in'h'c' sinh^(3:/2) 

^ ^''^ =7ir^r ^ Ia ^ ° if esc 6 - z esc 6 , B3 

where b is Galactic latitude. Integrating e(r) over the disk volume gives the total luminosity 
of the Milky Way 

L = 2H J 27irdr eoe^''/^° = ^TiHrleo. (B4) 

Combining Eqs.BS and B4, we can obtain the total luminosity of Milky Way from the 
observed Galactic emission 



L = ry/-> . . . I sin b\ [^(6; r) + ^^(-6; r) 



27r/i3c2 sinh2(a;/2) 
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(B5) 



Figure 7 in p gives the cosecant law for the Galactic emission in different WMAP bands. 
Using this information in Eq.(B5) (after conversion into thermodynamic units) gives the 
luminosity of the Milky Way in WMAP bands 

L*Q = 1.7 X 1037 ergs-\ 
L*y = 3.0 X 10^7 ergs~\ 
and = 1.0 x lO^^ ergs-^ (B6) 

To confirm these values, we can use Eq.(B2) and the observed integrated flux of the 
Andromeda (M31) galaxy in the WMAP maps to obtain its microwave luminosity 

Lm3i,q = 2.1 X 10^7 ergs-i, 
Lm3i,v = 5.3 X 10^7 ergs-\ 
and Lm3i,w = 1-6 x 10^^ ergs^^ (B7) 

We see that these values are larger, but within 50%, of the Milky Way microwave luminosi- 
ties. 
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